Pricing of barrier options by marginal functional quantization 

Abass SAGNA * 

■ Laboratoire d' Analyse et de probabilites de 

C^l , rUniversite d'Evry Val d'Essonne & 

q : ENSUE 

Q 
in 



Abstract 

This paper is devoted to the pricing of Barrier options by optimal quadratic quantization 
Oh I method. From a known useful representation of the premium of barrier options one deduces an 

■ algorithm similar to one used to estimate nonlinear filter using quadratic optimal functional quan- 

ta I tization. Some numerical tests are fulfilled in the Black-Scholes model and in a local volatility 

model and a comparison to the so called Brownian Bridge method is also done. 

1 Introduction 

> 

^ I Consider a fixed time horizon T, which will be typically the maturity of the option in a financial 

O ■ model, and let {Q,,J^,F) be a probability space (modeling the randomness of the market) with a 

filtration T = {Tt, 1 < t < T} satisfying the usual requirements. The probability P is supposed to 
^ I be the probability in the 'real world' in opposite to the risk neutral probability. 

^ ' Consider that the stock price process (Xt)te[o.T] satisfies the following time homogenous stochastic 

differential equation (SDE) 



dXt = h{Xt)dt + a{Xt)dWt, Xq = x E M, (1.1) 

where (Wf)t6[o,T] denotes a one-dimensional Brownian motion defined on the probability space 
P); b : M — M and o" : M — M are continuous functions satisfying the global Lipschitz 
and linear growth conditions: 

\bix) - b{y)\ + \aix) - a{y)\ < C\x - y\ (1.2) 

and 

\b{x)\ + \a{x)\<C{l + \x\), (1.3) 

for every t E [0,T] and for every .x, y E M. The filtration considered here is the natural filtration of 
the brownian motion completed by the P-null sets. 
It is known that under the above ass ump tions on the coefficients of the diffusion there exists a unique 



strong solution for the SDE (see e.g. Ill2l . ll6l1 ). The uniqueness of the solution is ensured by the global 
Lipschitz assumption_ (ll.2p whereas the linear growth assumption ()1.3p guaranties that this solution 
do not explode (see lllql for more details). 



*E-mail: abass.sagna@gmail.com. This research! is supported by the "Chaire Risque de Credit" of the French Banking 
Federation. 
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The first workable model for 'rational' market pricing of traded options have been proposed by 
Black-Scholes in 1973 and extended by Merton in the same year. In the Black-Scholes model the 
economics consists on two assets: the stock price with dynamics as the previous SDE with b{t, x) := 
fix and a{t, x) := ax, and a zero-coupon bound of constant interest rate r and maturity T. 

Moreover, we know that under arbitrage free and completeness assumptions, the discounted price 
at time t of any European contingent claim is uniquely determined and is the expectation, under a 
probability P called risk neutral probability, of its discounted payoff (a functional of the price process 
(^t)ie[o.T] which may depend on all the trajectory of the process), given all the information available 
up to time t. If Vt is the value of the option a time t and if h denotes the payoff at the maturity, then 

Vt = e-^(^-*)E(/i| J-i), 

where E is the expectation under P, so that the price at time is 

Vo = e-''^E(/i). (1.4) 

Our aim in this work is to estimate such an expectation for a class of path-dependent payoffs: barrier 
options, by optimal quantization method. We consider here a class of exotic options whose payoff 
depend on both the value of the underlying asset at the maturity and its maximum or its minimum 
over [0, T]. This means, payoffs h of the form 

h = F(XT, sup Xt) or h = F{XT, inf Xt). 

t6[0,T] te[0,T] 

When the payoff can be decomposed as 

h = (/j(Xr)l{sup,g[o,y] Xt&i} or (^(XT)l{inf,g[o.Ti Xt&i} 

where / is an unbounded interval of M, one speaks about barrier options. This last class is a particular 
case of payoffs of the form 

h = ip{XT)l{ro{X)>T}, 

where T£,{X) is the exit time of a domain D C M'' by a d-dimensional underlying asset X = 
{X\...,X<i). 
Here are some useful definitions. 

Definition. The option is said to be an up-and-out option if it knocks out when the price of its under- 
lying asset crosses a specified value. It is said a down-and-out option if it has barrier below the initial 
asset price and knocks out if the underlying asset price falls below the barrier. 

The payoff of an up-and-out call expiring at time T, with strike price K and up-and-out barrier L is 
given by : 

{Xt - K)^\i^ sup Xt<L} 

te[o.T\ 

and the payoff of a down-and-out call barrier option with maturity T, strike K and barrier L is given 
by 

(Xt - -fSr)"^l{ inf Xt>L}- 

t6[0.T] 

The payoff of put options are defined similarly with {K — Xt)'^ in place of {Xt — K)'^ . 
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Note that closed formulas are available for the price of such options in the Black-Scholes framework, 
see 10]. But this no longer holds when we move out from the Black-Scholes framework. So that 
we are led to estimate the prices by some numerical procedures. One of the used methods is the 
regular Brownian bridge method (see e.g. [1]). It provides approximation formulae of the price 
of barrier options using diffusion bridge methods. This leads to useful forms to approximate these 
prices from recursive formulas (already pointed out in ifij]) similar to an algorithm used in 1 18] to 



estimate nonlinear filter by optimal quantization method. One difference of our setting with respect 



to the one of IllSn is that our algorithm involves non-regular functions. Furthermore, if we consider 
local volatility model in the previous setting, one way of processing the algorithm is to use Lloyd's 
algorithm (or stochastic algorithms) to compute the optimal grids and the transition probabilities. But, 
because of the irregularity of functions appearing in our context, one must increase the grid sizes of 
the marginal quantization of the price process to obtain good approximations of the prices. It is clear 
that this will be very time consuming to use Lloyd's algorithm to compute grids sizes, and, this also 
depends to the parameters of the model. Moreover, the marginal quantized process is not a Markov 
chain and, for numerics, it is forced to satisfy the Markov property. 

In this work, we propose a procedure based on (quadratic) marginal functional quantization method. 
It consists first in considering the ordinary differential equation (ODE) resulting to the substitution of 
the Brownian motion appearing in the dynamics of the price process (11.11 ) by one quadratic quanti- 
zation of the Brownian motion. Then, constructing some "good" marginal quantization of the price 
process based on the solution of the previous ODE's, we show how to estimate the premium of barrier 
options from a recursive formula similar to an algorithm used to estimate nonUnear filter using optimal 
quantization method. Note that by construction, the marginal quantized discrete process is a Markov 
chain. Furthermore, because this procedure is based on the quantization of the Brownian motion, it 
does not depend on model parameters and price estimates are obtained in few seconds (at most in 6 
seconds and sometimes instantaneously, for considered examples). Numerical simulations are per- 
formed in the Black-Scholes model and in the local volatility model called a pseudo CEV model. A 
comparison with the Regular Brownian Bridge method show that the former method may some times 
be faster and competitive with respect to the last one. 

The paper is organized as follow. Since in a general setting, the estimation of the prices requires 
paths discretization of the process, we will recall in Section [2] the Euler scheme and some relevant 
convergence rate. Then, we will see in Section |3]how to derive the price estimates from (continuous) 
Euler scheme. The obtained formulas ai^e well known and are moreover in a useful form to apply an 
algorithm similar to that used in nonlinear filtering estimation via optimal quantization. The algorithm 
and the relevant error are given in Section |4] This algorithm involves the marginal quantization of 
the stock price process and, in Section [5] we show how to construct such a process from a basic 
construction of functional quantization of a diffusion process. We end by some numerical experiments 
where we compare our method with the regular Brownian bridge method in the Black-Scholes model 
and in the pseudo CEV model. 

2 Euler Scheme 

Consider a one-dimensional Brownian diffusion process (^t)jg[o,T] > solution of the following stochas- 
tic differential equation 

dXt = b{Xt)dt + a{Xt)dWt, = x G M (2.1) 
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where 6 : M — )• M, cj : M — )• M are continuous functions satisfying conditions (11.21 ). (11.31 ) and 
(^t)te[o,T] denotes a one-dimensional Brownian motion defined on ($7, T ^ P). 

Let us divide the set [0, T] into n subsets of length T/n and set for every A; = 0, . . . , n, = 
The stepwise constant Euler scheme is defined by 



Xt,^, =Xt, +6(XtJ - + Xo = x, A: = 0,...,n-1 (2.2) 

where {Zk)i<k<n is a sequence of i.i.d random variables distributed as J\f{0; 1). 

For every t G [0, T], set i = if t G [t^, tfe+i), A; = 0, . . . , n — 1. A natural extension of the 
discrete Euler scheme is the continuous Euler scheme defined for every t G [0, T] by 

Xt=Xt + b{Xt) {t-t) + a{Xt) {Wt -Wt), Xo = x 

which satisfies the SDE 



Xt=x + [ b{Xs)ds+ [ a{Xs)dWs. 
Jo Jo 



The above paths discretization methods generate some errors which estimates are given in the follow- 
ing results (see e.g []). 

> Strong error rate. Assume b and cr satisfy for every a G (0, 1), 

VtG [0,r],Vy,zGM^ \b{s,y) - b{t, z)\ < C{\t - + \y - z\). (2.3) 

Then, 

(a) for every p > 0, for every n > 1, 



II sup \Xt-Xt\\\p<Ct,a,pe'^^''-'^{l + \x\) 
ie[o,T] V"- 

(6) for every p > 0, for every n > 1, 

sup \Xt - Xt\\\p < Cb^^^p e''-''^'-'P{l + \x\) 



te[o,T] 



n 



> Weak error. We recall some weak error estimates for path-dependent options (we refer e.g. ifioll for 
the proofs). Let 

D([0,r],]R'^) := |e : [0, T] ^W^, cadlag} . 
If F : ][D([0, T], M'*) M is a Lipschitz functional for the sup norm, that is, 

\F{i)-F{i')\<CF sup \m-at)\ 

t6[0,T] 



then 



and 



m{{Xt)te[o,T\) - EF((XOt6[o,T])| < ^ (2.4) 



' loo;' 71 

\¥.F{{Xt)te[o,T])-^F{{Xt)te[o,T])\ < C\l^. (2.5) 



n 
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On the other hand, if a domain D has a smooth enough boundary, 5, cj G C^(M) and a uniformly 
elhptic on D : 3ao > 0, Vx G M a^{x) > ctq, then, for every bounded measurable function / 
satisfying d(supp(/), dD) > 2e > 0, 

and 

where n is the number of discretization steps and t{Y) is the exit time of the process Y from the open 
set D, i.e 

t{Y) = mi{t G [0,T], Yt G D^}. 

Then the convergence rate is of order for the continuous Euler scheme and of order n"^/^ for the 
discrete one. 

3 Approximation of knock out option prices using diffusion bridge 

According to the convergence rate for the continuous Euler scheme we would like to estimate the price 
of path-dependent options by replacing the asset price process (^i)te[o,T] by its continuous Euler 
process (^)t(=[o,T] (ll-4b . Then, given values of the process (Xt) at discrete times tk, k = 0, . . . , n, 
one deduces formulas integrating the useful information which is the probability that the barrier is not 
knocked over the time interval [0, T]. Remark that this information is lost when replacing X by the 
discrete Euler process X in (11.41) because we do not known if whether or not the barrier is knocked 
between time intervals (tfc, t^+i), k = 0, . . . , re — 1. But, integrating this information in the former 
case requires the knowledge of the distributions of the maximum and the minimum of the continuous 
Euler process {Xt) over the time interval [0, T], given its values at the discrete time observations t^. 

Proposition 3.1. We have 

C{maxXt\Xt,^=Xk,k = 0,...,n) = C{ max G'^ (Uk)) (3.1) 

and 

^^m^] =x,,k = 0,...,n)= /:(^^min _^ F.-',.,^AUk)) (3-2) 

where {Uk)k=o,...,n-i f^f^ i-i.d random variables uniformly distributed over the unit interval, 
and F~y are the inverse functions of the conditional distribution functions Gx,y and Fx^y defined by 

1 _ e TaH^) j l|„>max(x,y)} 

and 

1 _ e T.'(^) \ l|„<n,in(x,s,)} • 

This result is proved using the independence property of the processes {Xt)t^[t^.^tk+i]^ ^'^^ ^ ~ 
0, . . . , re — 1, given the Xt,, = Xk, and the knowledge of the distribution of the supremum (and the 
infimum) of brownian bridge diffusion over time intervals (tfc, t^+i), with end points x^ and x^+i- 

From the above proposition we deduce general formulas making a connexion between the expecta- 
tion of a functional of both the terminal value Xt of the process (Xt) and its maximum (or the min- 
imum) over the time interval [0, T]. From now on we make the abuse of notation := Xt^, \/k G 
{0,...,re}. 
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Proposition 3.2. (a) Let f be a real-valued non negative function defined on such that f{x, •) is 
a nonnegative function satisfying 



sup E/(2;, max Xt) < +00. (3.3) 

x>o te[o,T] 



Then 



E/(Xt, max Xt) = E/(Xt, 0) + E / 1 " TT 1 x.i^^ d-JiXr, z). (3.4) 
te[o,T] 7o \ ^^^^ / 



Likewise if 



supE/(3;, mill Xt) < +00) 
x>o te[o,T] 



/■+00 / " \ _ 

E/(XT,^min X^VO) =E/(Xt,0)+E fj] (l - F;f^_^^jf^.(^)) 1 '^./(^T, (3.5) 



(6) If furthermore fooix) ■= lim f{x,y) < +00 for every x > 0. Then 

3/— >+oo 



EfiXr, max X^) = E/oo(Xt) - E / TT 4/(Xt, (3.6) 



and 



r+00 / " \ _ 

E/(XT,^min V 0) = E/oo(Xt) - E ( 1 - ^ (l " 1 ^^./(^T,^). (3.7) 

This proposition follows from the following lemma. 

Lemma 3.1. (a) Let Z be a positive random variable and let g be a nonnegative function with finite 
variation (on compact sets) such that 



Then 



'(0,+oo) 

(6) If furthermore Qoo ■= lim g{x) < +00 then 



E I / \dg\ < +00. (3.8) 
]o,z] / 



Eg{Z) = g{0) + [ P(Z > z) dg{z). (3.9) 



Eg{Z) =goo- I z) dg{z). (3.10) 

J{0,+oo) 
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Proof, (a) We have 



It follows that 



g{Z) = 5(0) + / dg{u). 

]0,Z] 



¥.g{Z) = g(0)+E / dg{u) 



= g{0) + F{Z> z) dg{z), 

J(0,+oo) 

the last inequality coming from Fubini's theorem; which can be applied owing to assumption (|3.8p . 
(b) Just use the fact that P(Z > u) = I - F{Z < u). □ 

Now we are in position to prove Proposition 13.21 

Proof of Proposition 13.21 (6) One deduces from ()3.ip that 

E(f{xn, max Xt)\Xk = Xk,k = 0,. . . ,n) =E(f{xn, max G^^., , , (t/fc))) 

where and the C/^ are defined like in (|3.ip . Then, applying Lemma ITT] (6) to the function 

g{z) = fixn,z) gives 

r+oo 

E{f{xn max G-^^^^^iUk))) = foo{xn)- F{ max C'l^^^^iUk) < z)d J {x^, z) 



foo {^r, 



r+oo /n-l \ 



vA;=0 



Xk,Xk+iiz)) dzf{Xn,z). 
\fc=0 / 



Consequently 



EfiXr, max Xt) = E(E(/(Xt, max Xt)\Xk = Xk, k = 0, . . . ,n)) 



r+oo / n \ 

E/oo(Xt) -E [JlGx,.,,x,iz)j dJiXT,z). 



The formula relative to the minimum is proved likewise by using (|3.2p in place of ()3.ip . 

(a) is proved like (6) by using Lemma ITT] (a) instead of Lemma 13.11 (6). □ 

Proposition l3.2l allows us to rewrite the estimates of the premiums of some usual exotic options (in 
particular barrier options) in a useful form in view of the optimal quantization approximation method 
as well as of Monte Carlo simulation methods. Let us mention that the following representations of 
the price estimates of Barrier options are well known even if the computational method used here to 
derive them is a little different. 
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Proposition 3.3. Let f{x) = {x — K)~^ and g{x) = {K — x)'^. 

(a) The price of an up-and-out put option expiring at time T with strike K and up-and-out barrier L 
is estimated by 

Pub := e-^^E((K - sup x,<L}) = e'^^K \9{Xt) f{ G^, ) • O-H) 

te[o,T] \ ^.^1 / 

(6) The price of an up-and-out call option expiring at time T with strike K and up-and-out barrier L 
can be approximated by 

Cub := e-'-^E((Xr - K)+\^ ^^^^^ = e^^^E ( /(X^) J] ) • (3-12) 

te[o.T] \ fc^-^ J 

(c) The price of an down-and-out put option expiring at time T with strike K and down-and-out 
barrier L can be approximated by 

PoB:=e-^^E((/^-XT)+l| x,>L})=e-'-^E(5(XT)n(l-F^,_,,x,(^))V (3-13) 

(d) The price of an down-and-out call option expiring at time T with strike K and up-and-out barrier 
L is approximated by tite following formula : 

Cob := e"^^E((XT - K)+l^ = e-'^^E f /(X^) J] (l - F^^_^,^^(L)) ) . (3.14) 

Note that the right hand side of Equations p.l ip . ()3.12p . p.l3p . p.l4p are obtained by re-conditioning 
Then, it follows from Jensen inequality that the corresponding variances are smaller than the variances 
induced by the left hand side of the same equations. 

Proof, (a) Let f{x) = {x — g{z) = 1{z<l} and set h{x, z) = f{x)g{z). Then itfoUows from 
ra that 



E((Xt - x,<L}) = ^hUXr) - E / f{XT) J] ^x^.^x,!^) dg{z). 

te[o,T] Jo \ fc^^ J 

Now Vx > 0, hoo{x) = and dg{z) = —Sl{z). Then 

Cub = e-^^E ( /(Xy) fl ) . 



fc=i 



The items (6), (c), (c?) are proved in the same way as (a). □ 

In the next section, we show how to estimate the previous prices by optimal quantization. We will 
give first an approximating algorithm and then, the induced en^or. Since this algorithm involves the 
quantization of the price process and its transition probabilities, we will point out how to construct a 
functional quantization of the price process and how to estimate its transition probabilities. 
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4 Estimation of the prices by marginal quantization 



The aim of this section is to propose an algorithm based on optimal quantization to compute the path 
dependent options pointed out in Section[3l To this end we will approximate by optimal quantizations 
some expressions of the form 

V:=E (^f{Xn)fl9k{Xk^i,Xk)j (4.1) 

where / is a bounded measurable function on taking values on M and gk{-,-) a measurable function 
on R'^ X which may depend on some real parameters like for Barrier options where it depends also 
on the barrier. 

4.1 The algorithm 



We will mainly refer to jlSn . where numerical solving of nonlinear filtering with discrete-time obser- 



vation have been performed by optimal quantization methods. The only change is that in our setting 



we will drop the dependance on the noisy observations (i.e the Y^'s following the notations in jlSy) 
because our problem of interest here is not a filtering problem. 
We define for any k = 1, . . . ,n, the bounded transition kernel Hk by 

Hkf{x) = E {f{Xk)gk{x,Xk)\Xk-i = x)= I f{y)gk{x,y)Pk{x,dy) (4.2) 

where Pk{x, •) = C{Xk = -l^fc-i = x). For convenience, we set 

Ho fix) = E(/(Xo)) = J f{x)fi{dx). (4.3) 

Now for any k £ {1, . . . ,n} set 

k 

TTkf = E{f{Xk)llgk{Xi-uXi)). 

1=1 

We have 

T^kf = E ^E(/(Xfe)J]5,(X,„i,X,)|Ji,_jj 

= E ^E(/(Xfc)5fc(Xfe_i,Xfc)|Xfc_i) J]5,(X,_i,X,)j 
= E ( Hk{f{Xk_,))l[g,{Xi_i,Xi)] 



It follows that 

TTfc/ = TTk-iHkf, k = l,...,n (4.4) 

so that 

V = 7rnf = {HooHio...o Hn)f. (4.5) 
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Then, to estimate V we need to approximate 7r.„. At this step, suppose that we have access to 
the quantization {X)t^ of the price process over the time steps tk,k = 0, . . . , n on grids Tk = 
{xl., . . . , x'^''} of sizes Nk, k = 0, . . . ,n (see further on for facts about quantization). 

Owing to equation ()4.4p our aim is to estimate the price using an approximation of the probability 
transition Pk{xk,dxk+i) of Xk+i given Xk- These probabihty transitions are approximated by the 
probability transition matrix := {p'^) of X^j^i given X]^: 

= P(Xfe = = 4_i); i = 1, . . . , iVfc_i, i = 1, . . . , Nk. (4.6) 

Then, following Equation (|4.2p . we estimate the transition kernel matrix by the quantized tran- 
sition kernel given by 

Hk = ^Hl^S,.^_^, k = l,...,n 
i=i 

where 

Hi' = gki4-i,4)pt ^ = 1, • • • , Nk-i, j = l,...,Nk. (4.7) 
For A; = 0, we set (owing to ()4.3p and to the fact that Xq = xq is not random) 

Ho = SxQ- 

We finally approximate 7r„ by 

nn = HooHio...oH^; (4.8) 
which in turn can be computed by the forward induction 

vro = ^0, TTfc = 7ffc_i#fc, k = l,...,n. (4.9) 

It follows that the price V = TTnf niay be estimated by summery 

V := ^nf. 

From the previous approach, we deduce the following estimations for options of interest using 
optimal functional quantization method. Set in this scope f{x) := {x — K)~^ and g{x) := {K — x)~^. 

> Up-and-out options. According to the forgoing we estimate the price of an up-and-out put option 

by 

Pub ■= e ''^vr„5 
and the price of up-and-out call option is approximated by 

Cub ■= e~^^7f„/ 

where ??„ is defined as in (|4.8p with the associated transition kernel 

Hl^=G^, ^j{L)p'i, i = l,...,iVfc_i; j = l,...,iVfc. 

> Down-and-out options. The down-and-out put option's price is estimated by 

PoB := e~'''^7fng 
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and the price of down-and-out banier call option is estimated by 

Cob ■= e'^^TTnf 

where for both cases 7f.„ is defined as in ()4.8p with the associated transition kernel 

= F , (L)p^^ i = 1, . . . , iVfc_i; j = l,...,Nk. 

Remark 4.1. One numerical advantage of this algorithm is that does not depend on the function 
/ appearing in (I4.5I ). Then, once 7f„ is computed we deduce both the call and the put price approxi- 
mations. On the other hand, considering Equation (14.71 ) one notices that as soon as X^-i reaches the 
barrier (for example, for the up-and-out option: there exists zq such that > L), then, H^f^ = for 
every i > io- For numerical computation, we may take account of this fact to reduce the computation 
time. 

4.2 Error analysis 

In order to have some upper bound of the quantization error estimate of vr/ we need the following 
assumptions (Al) and (A2) : 

(Al) The transition operator Pk{x, dy) of given = 1, . . . , n are Lipschitz. 

Recall that a probability transition P on is C-Lipschitz (with C > 0) if for any Lipschitz function 
/ on with ratio [f]Lip, Pf is Lipschitz with ratio [Pf]Lip < C[f]Lip- Then, one may define the 
Lipschitz ratio [-PJljp by 

[P]Lip = sup {K-pr—^, f a nonzero Lipschitz function } < +oo. 

[/ J Lip 

Then if the transition operators Pk{x, dy), A; = 1, . . . , n are Lipschitz, it follows that 

[P]Lip ■■= max [Pk]Lip < +00. 

fc=l,...,n 

(A2) It consists on the following two assumptions. 

(i) For every k = 1, . . . ,n, the functions gk{-, •) are bounded on R*^ x and we set 

Kg := max ||5fc||oo 

k=l,...,n 

(a) For every k = I, . . . ,n, there exist two constants [gl]Lip and [gl]Lip so that for every x, x', x, x' e 

\gk{x,x') - gk{x,x)\ < [gllup \x-x\ + blliip W - x'\. 

Theorem 4.1. Under Assumptions (Al) and (A2) we have for every bounded Lipschitz continuous 
function f on and for every p > 1, 

n 

\7Tnf - ^nf\ < C^(/,p) \\Xk " Xk\\p (4.10) 
fc=0 

with 

Ct{f,p) = (2 - 62,p) [ukhip + K^-' \\f\U[gUl]L^p + [gUi]up). 
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Proof. The proof follows from the proof of Theorem 3.1 in [ Ts^ by dropping the dependency on the 
noisy observations (yi, . . . , yn) following the notations of the authors. □ 

Now, let us come back to the construction of the quantized price process (^fc)fc=o,...,n- We show in 
the next section how to construct this process after making a short background on product functional 
quantization of gaussian processes, in particular, of brownian motion. 



5 Marginal functional quantization of the price process 

Before dealing with the construction of the marginal functional quantization of the price process, we 
make some background on functional quantization of gaussian processes. 



5.1 A brief overview on functional product quantization of gaussian processes 

We remind first some basic notions about optimal vector quantization. It is a process of approximating 
a continuous range of values or a very large set of discrete values by a relatively small set of discrete 
values. Rigorously speaking, the L''-optimal quantization problem at level n for a R'^-valued random 
vector X lying in U{Q., P) consists in finding the best approximation of X by a Borel function of 
X taking at most n values. This problem can be reads as 

en,r{^) = inf{||X -X"||^,a C M"^, card(a) < n} 

inf (I d{x,aYdP{x) \ ' . (5.1) 

aCM"^ \JM.d ) 
card(ci)<n 

where X" = X^^jg^ Q^l{X6Ca(a)} quantization of X on the grid a and (Ca(a))aeQ corresponds 
to a Voronoi tessellation of (with respect to a norm | • | on M"^), that is, a Borel partition of M 
satisfying for every a G a, 

Ca{oL) C {x G E'^ : |x - a| = min \x - 

The quantity en^ri^X^ is called the L'^-mean quantization eiTor. This eiTor decreases to zero at a 
j^-i/'^.j-ate as the size n of the codebook a goes to infinity. This convergence rate has been investi- 
gated in ^\ and lis!] for absolutely continuous probability measures under the quadratic norm on 
and studied in great details in |jl Ij] under an arbitrary norm on for absolutely continuous measures 
and some singular measures. Very recently, optimal vector quantization has become a promising tool 
in Numerical Probability owing to its ability to approximate either expectations or more significantly 
conditional expectations from some cubature formulas. This faculty to approximate conditional ex- 
pectations is the crucial property used to solve some problems ernerging in finance as optimal stopping 
problems (pricing and hedging American style options, see 10, 120I1 . stochastic control problems (see 
BE3) for portfolio management, nonlinear filtering problems (see |18, 21] and |3l for an application 
to credit risk). 



A rigorous extension of optimal vector quantization to functional quantization is done in \\% where 
the vector quantization problem is transposed to random variables taking values in an infinite dimen- 
sional Hilbert space, in particular, to stochastic processes (Xt)(g[o 1] viewed as random variables with 
values in i^^([0, 1], dt). Many others works have been done in this direction as e.g. From the 



numerical point of view, it is pointed out in ll20n how a Gaussian process can be quantized using 



12 



Karhunen-Loeve product quantization based on product quantization of Gaussian random variables 
coming from the Karhunen-Loeve expansion of the given Gaussian process. A closed formula for the 
distribution of the quantization of the Gaussian process (in particular for Brownian motion) is derived 
and some applications has been successfully performed in Finance, namely, in the pricing of vanilla 
and Asian call options in Heston model. 

To recall some basic results about functional quantization suppose that {H, {•\-)h) is a separable 
Hilbert space and let X : {Q.,A,¥) H he square integrable if- valued random vector with distri- 
bution Fx defined on {H, Bor{H)) where Bor{H) stands for the Borel ci-field. Let || • || denotes the 
L|^(f],P)-norm defined by \\X\\l = E{\X\]j). 

Let X := {xi, . . . ,Xn} G be an ?i-quantizer and let X^ be the quantization of X on the grid 
X defined previously, where the Voronoi tessellation (Ci(x))i<i<„ induced by x satisfies for every 
i G {1,... ,n}, 

Ci{x) C {y £ H : \xi - y\H = min \xj - yln}- 

l<j<n 

The quadratic quantization problem consists of finding an optimal quantizer x G (if any), means, 
an n-quantizer which minimizes the quantization error \\X — X^\\2 over H"^. From the numerical in- 
tegration viewpoint, finding an optimal quantization may be a difficult problem and we are sometimes 
let to find some 'good' quantization X^ which is close to X in distribution, so that for every Borel 
function : if i— M, we can approximate KF{X) by 

n 

EF(X^-) = J2Fixi)Fxia{x)). (5.2) 

Then if we have access to both the n-quantizer x = {xi, . . . ,Xn} and the distribution associated 
to X^, (Px(Cj(x)))-|^^^^^^, the estimation of ¥,F{X) using Equation (15.21 ) is straightforward. The 
induced error depends on the regularity of the functional F and here is some error bounds. Suppose 
that X e L^(0, P) and let F be a Borel functional defined on H. 

1. If F is Lipschitz continuous with Lipschitz constant [F]Lip then for every n-quantizer x, 

\EF{X) - EF(1"')| < [F]up\\X - l^lla 
so that if {xn)n>i is a sequence of quantizers satisfying lim \\X — X^"\\2 = 0, then X^" 

~ n— >-oo 

converge in distribution to X. 



2. If F is differentiable on H with an ^-Holder differential DF, 6 G (0, 1], then for every optimal 
n-quantizer x, 



\EF{X) - EF(X^)| < \DF\e\\X - 



Some others bound are available (we refer to 112011 for more detail). Now let us say how to get some 
'good' quantizers for Gaussian processes to make sense the previous errors bounds. We consider here 
a centered one-dimensional L\ := L^([0, T],dt)-vahied Gaussian process X satisfying 



E 



X\l2 = [ E{X'^)ds < +00. 
^ Jo 



The process X admits the following representation in the Karhunen-Loeve basis (see e.g. 111411 ) 



X{u) ^ Yl V^kU^^k - a.s. (5.3) 



fc>i 
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where the sequence {Ck)k>i defined for every > 1 by 

^Var((X|eJ)) 

is a sequence of i.i.d J\f{0; 1) -distributed random variables. Owing to the expansion (I5.4I ). a natural 
way to produce a functional product quantization of a Gaussian process in of size at most is to 
use a product quantizer of the form 

^^^ = i;v^et"^'ef(t) (5.5) 

k=l 

where is an optimal A'^fc -quantization of and cIn := NiX - ■ - x Nl < N, with A'^i, . . . , Nl > 

2. An quadratic optimal product A^-quantizers also noted X^"''^'' is obtained by solving the optimiza- 
tion problem (see (l^ for more detail): 

min{||A:-X('^'^)||2, dN = Nl X ■■■ X Nl < N;Ni,...,Nl>2;L>1}. (5.6) 

We suppose from now on that the previous optimization problem can be solved, at least numerically, 
and that the optimal L-tuple still be denoted by A'^i, . . . ,Nl. Then, numerical computation of a 
Gaussian process X is possible as soon as we have numerical access to the eigensystem {e^ , \n), 
which, for the Brownian motion (VFf)fg[o.r]> has a closed formula: 

er(t):=y|si„(,(t- 1/2)1) X,:={-^^y.t>l^ 

So, the one-dimensional quadratic optimal product quantizer , at level N, of the Brownian motion 

(^i)t6[o,T]> is defined by 

<,.... At) = E ("^^^ - 1 < < 1 < /c < L, 

where x^^*) = {x^'', . . . jX^^} is the optimal quantization of the J\f{0; 1) of size N^ and = 
Y[k=i is an optimal integer solving Problem ( 15.61 ) (with respect to the Brownian motion). Remark 
that for every t G [0,T], the marginal quantizer ^^(t) is of size d^- For numerics, a whole 
package of product A^-quantizers of the standard Brownian motion are available at www . quant i ze . 
maths- fi . com. We move now to the construction of the quantized price process. 

5.2 Marginal functional quantization and transition probabilities 

Recall that the continuous Euler price process evolves following the SDE 

dXt = b{Xt)dt + a{Xt)dWt, Xo = X. (5.7) 

Let {a^)N>i, with for every A^ > 1, {t) = {ai {t), . . . ,a^(t)}, be a sequence of optimal 
product A^-quantizers of the Brownian motion and let 

dN 

t?/'^-) = a^(t)l{H/,6C™(a-(t))}, t e [0,T], 

m=l 
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be the marginal functional quantization of the Brownian motion. It is known that the sequence 
{a^)N>i is rate-optimal, means, 

||Vr-T?"'^||2 = 0((logiV)-^/2). 
Consider the sequence = {x^)m=i,...,dff, N > 1 of solutions of the ODE's 

x^{t)=x+ f [b{xZ{s))-laa'{xZ{s))]ds+ f a{xZ{s))daZ{s), m = l,...,dN (5.8) 



^ Jo 

t = tfc if t G [tk,tk+i), k = 0, . . . ,n — 1, and define the marginal functional quantization of Xt over 
the grid x^(t) = {xf(t),...,x^^(t)}by 

^t^ = ^m(Ol{XteC„,(x^(t))}i (5-9) 
m=l 

SO that Equation ( 15.81 ) can be written as 

Xi^ = x+ f [b{X^)-\aa'{X^)]ds+ j\{X^)dWf-\ 
Jo ^ Jo 

Recall that the process (Xt^ ) is a Mai^kov chain. Then, since by construction o(Xt^ , k = 0, . . . ,n) = 
a{X^, k = 0, . . . , n), the discrete process {X^)k=o,....n is a Markov chain. On the other hand, 
since the additional term ^cra' appears in the ODE (this correction term can be dropped by consid- 
ering the stochastic integral in (15.71 ) in the sense of Stratonovich integral for ^) convergence 



investigation tools, see 11141] '). we must make the supplementary assumption that a is continuously 
differentiable with bounded derivative to guaranty the existence and the uniqueness of the solution. 
Now, given the quantization process {X'^), to complete the estimation of the price of banier options 
following the introduced algorithm in Section 14.11 it suffice to compute the transition probabilities 
appearing in (14.61 ). Notice that all our grids x^{tk) are of size d^ = Ni x . . . N^. To define cor- 
rectly the Voronoi cell associated to the grids x^{tk) we will consider that for every time step tk, 
x^{tk) = {xi{tk), . . . ,x^^{tk)} is a descendent ordered set. The computation of the transition 
probabilities will be made differently according to the following two situations. 

> TIte cumulative distribution function F{-;x) of the conditional law of Xt given Xs = x is known 
for every s < t. For example, this is the case in the Black Scholes model where F{-;x) is the 
cumulative distribution function of the lognormal distribution. In this case, since following Equation 



( 15.91 ). cr{Xt^, , k = 0, . . . ,7i) = (t(Xj^, k = 0, . . . ,n), the probabilities are estimated by 

« i^(xf+(tfe);xf (ifc_i)) - F(xf_(tfc);xf (tfc_i)), (5.10) 
with for every /c = 0, . . . , n — 1, 

Xj+{tk) •■= ; := 2^ ; j = l, . . . ,c?Ar - l; 



In fact, we have for every A; = 0, . . . , n — 1, 



P{Xk+i G Cjix^{tk))\Xk^i G Ci{x^itk-i))) 
P{Xk < xf^{tk)\Xk^i G Ci{x''{tk^i))) -P{Xk < xf„{tk)\Xk^i G C,(x^(tfc_i))). 
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Afterward, we have for every z > 0, 

and considering the numerator in the right hand side of the previous equation we have 

P{Xk < z;Xk-i G Ci(x^(tfe_i)))= ( / P{Xk G dx\Xk-i = y)dPx,_Ay)] dx 

= [ F{z-y)dPx,_,{y) (5.11) 

^F{z-xf{tk-i))P{Xk^i e 



The last quantity is the approximation of the right hand side of (15.111 ) by optimal quantization with one 
grid's point, considering that {x^^ (tk-i)} is the quantizer of size one of the random variable X^^.^ 
over the Voronoi cell Cj((x^(tfc„i)). 

I> The cumulative distribution function F{-,x) of the conditional law of Xt given Xg = x is unknown. 
In this case, considering the (discrete) Euler Scheme of the price process (see (12.21 )) we estimate F by 
the cumulative distribution function F of the M {rri}^; af,) with 

^ ^ T ^ T 

ruk = Xk-i + hitk,Xk-i) — ; al = a'^ (tk, Xk^i) —, 

n n 

so that for every i,j = l,..., d^, 

p'l « F(xf+(tfc);xf (tfc_i)) - F{x^_{tk);xf{tk-i)) (5.12) 

where the xj^(t/^.) and x^_{tk) are defined as previously. 

Notice that since the eiTor bound of the filter estimate in ( 14.101 ) involves the marginal quantization 
error: \\Xt^. — X^^\\2, one must deduce this error from the above construction. We know that the 
sequence of non- Voronoi quantization (X^ )n>i defined for every > 1 by 

m=l 



is rate-optimal in L% (17, P) for p G [1,2): \\\X - X"" |^2 |L = Oi {log N)-^l'^) (see UJ]). One 
theoretical challenge will be to compute the convergence rate for the mai^ginal functional quantization 
error. 



6 Numerical illustration 

We deal with numerical experiments by considering an Up-and-out call option in the Black-Scholes 
model and a local volatility model already considered in jlSll and called pseudo CEV model. Recall 
that in the Black-Scholes framework the stock price process {Xt) is modeled by the following SDE 
(under the risk neutral probabihty P) 

dXt = rXtdt + aXtdWt, Xq = xq (6.1) 
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where r is the interest rate, a the volatility and W a brownian motion under P. For the pseudo CEV 
model, the dynamics of the stock price process is ruled by the following SDE (under the risk neutral 
probability) 

dXt = rXtdt + ^Xf^^L=dWt, Xo = XQ (6.2) 
V 1 + 

for some 6 € (0, 1) and G (O,]?],]? > 0. The parameter r still be the interest rate and (t{x) := 
corresponds to the local volatility function. We notice that for a fixed 6 G (0, 1), if the initial 
value of the stock process Xq is large enough then the pseudo CEV model is very close to the CEV 
model 

dX't = X[{rdt + ^{X[f-^dWt). 

In particular, for numerical tests we will consider that -d « (tXq~^ where a denotes the regular 
volatility. The only "aim of the really" rough calibration is just to deal with reasonable values to 
obtain prices close to those given by the Black-Scholes model. For all the experiments we set the 
interest rate r equal to 0.15. The maturity is set to T = 1, the initial value of the stock process 
Xq = 100 and 5 = 0.5 (in the local volatility model). For numerics, the solution of the ODE given 
in (15.81 ) is approximated by a sixth order Runge-Kutta scheme and marginal quantizations are of size 
dfyf = 966 (corresponding to the optimal decomposition A^i = 23, N2 = 7, = 3, = 2, for the 
problem ^M, see Q). 

In the Black-Scholes model we compare the prices computed from the quantization of the contin- 
uous Euler process using (15.81) and ( 15.91 ) (which prices are referred by QEP prices) with the regular 
Brownian Bridge method (RBB prices), given the true prices obtained from a semi-closed formula 
available in |@]. The regular Brownian Bridge (RBB) method is some efficient method to compute 
expressions like 

E/(Xt, sup Xt) or E/(Xt, inf Xt), 
te[o,T] *e[o,T] 

based on Proposition 13. H and consisting (for example for the estimation of E/(Xt, sup^g[o,T] ^t)) 
the following steps : 

Set Sf = 0. 

for m = 1 to M 



Simulate a path of the discrete time Euler scheme (X^™^) and set Xk = xl:"^\ k = 



,n. 



Simulate T^™) := max (G^,,^,^ J-i([/f '^), where (f/^ )i<fc<n areiid with^/([0, 1]) -distribution. 

compute /(X^"'\r(™)). 

Compute SL := /(xi!"\r(-)) + Si_^. 



end. (m) 
with 



{G,,y)-\l-u) = -{x + y+y/{x- 2/)2 - 2Ta^x) log(n)/n) , uG (0,1). 
Then for large enough M, 

5/ 

EfiXT, sup Xt) « -f . 
te[o,T] J-VJ- 
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The numerical results are depicted in Table |2] and Table [3] for varying values of the barrier L and the 
volatility a. The number M of Monte Carlo simulations is set to 10^. For the quantization methods, 
the computation times (QEP c.t.) varies from to 3 seconds when n = 10 and, from 1 to 6 seconds 
when n = 20, increasing with the barrier as pointed out in Remark 14771 However, for the RBB, 
the computation time is of 2 seconds when n = 10 and of 5 seconds when n = 20. The obtained 
results show that the quantization method may sometimes be competitive with respect to the regular 
Brownian Bridge method, specially for small number of time discretization steps n. 



a = 0.07, n = 10 


L 


True prices 


RBB prices 


RBB vai-. 


QEP prices 


QEP c. t. 


105 


0.034 


0.035 


0.086 


0.035 


<ls 


110 


00.59 


00.60 


2.942 


00.59 


Is 


115 


02.58 


02.62 


15.80 


02.59 


2s 


120 


06.01 


06.11 


33.54 


06.03 


2s 


125 


09.58 


09.89 


41.76 


09.60 


2s 


130 


12.07 


12.13 


43.09 


12.08 


3s 



Table 1 : Up-and-out Call prices from the quantization method (QEP prices) and the regular brownian bridge 
method (RBB prices) in the Black-Scholes model for r = 0.15, cr = 0.07, 966, T = 1, K = 100, 

Xq — 100, n = 10 and for varying values of the barrier L. QEP c.t. is the quantization method computation 
time. 



cr = 0.07, n = 20 


L 


True prices 


RBB prices 


RBB vai-. 


QEP prices 


QEP c. t. 


105 


0.034 


0.035 


0.086 


0.034 


Is 


110 


00.59 


00.60 


2.942 


00.59 


3s 


115 


02.58 


02.59 


15.80 


02.59 


4s 


120 


06.01 


06.05 


33.54 


06.02 


4s 


125 


09.58 


09.64 


41.76 


09.59 


5s 


130 


12.07 


12.10 


43.09 


12.08 


6s 



Table 2: Up-and-out Call prices from the quantization method (QEP prices) and the regular brownian bridge 
method (RBB prices) in the Black-Scholes model for r = 0.15, a = 0.07, dN = 966, T ^ 1, K = 100, 
Xq = 100, n = 20 and for varying values of the barrier L. 

For the local volatihty model we compare the QEP prices with the prices obtained from regular 
Brownian bridge method. Numerical results are depicted in Tables |4] and [5] for different volatilities 
and for different values of the banier. Our reference prices are computed from the regular Brownian 
bridge method with 10'' Monte Carlo simulations and 100 times discretization steps. We remark that 
the quantization method is competitive (because it is faster with the same precision) with respect to 
the regular Brownian bridge method when the barrier is closed to Xq. But, the QEP prices become 
less precise when the volatility and the barrier increase (for example when cr = 1.0 and L = 130, 
where the absolute enor is of 2% with respect to RBB method). This might be due to the additional 
eiTor coming from the estimation of the transition probabilities by formula (15.121 ) since the conditional 
law is not known. 

Notice that, since we have used sixth order Runge-Kutta scheme to approximate solutions of the 
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a = 0.1, n = 20 



L 


True prices 


RBB prices 


RBB var. 


QEP prices 


105 


0.029 


0.029 


0.067 


0.029 


110 


00.42 


00.43 


1.933 


00.42 


115 


01.70 


01.72 


10.46 


01.71 


120 


03.95 


03.98 


26.42 


03.97 


125 


06.70 


06.76 


43.82 


06.72 


130 


09.31 


09.38 


57.19 


09.34 



Table 3: Up-and-out Call prices from the quantization method (QEP prices) and the regular brownian bridge 
method (RBB prices) in the Black-Scholes model for r = 0.15, a = 0.1, cIn = 966, T = 1, K = 100, 
Xq — 100, n — 20, and for varying values of the barrier L. 



1} = 0.7, n = 20 



L 


Ref. Price 


RBB price 


RBB var. 


QEP price 


QEP c. t. 


105 


0.034 


0.034 


0.085 


0.034 


2s 


106 


0.074 


0.074 


0.222 


0.074 


2s 


107 


00.14 


00.14 


0.496 


00.14 


2s 


110 


00.59 


00.59 


2.949 


00.59 


3s 


111 


00.86 


00.86 


4.636 


00.86 


4s 


112 


01.20 


01.20 


6.841 


01.20 


4s 


115 


02.64 


02.66 


16.34 


02.66 


5s 


120 


06.25 


06.29 


34.47 


06.30 


6s 


125 


09.98 


10.00 


41.73 


10.02 


7s 


130 


12.44 


12.45 


41.65 


12.46 


8s 



Table 4: Up-and-out Call prices from RBB and QEP methods in the local volatility model. Model parameters; 
r = 0.15, (5 = 0.5, ■& = 0.07, = 966, T ^ 1, K ^ 100, Xo = 100 and for varying values of the barrier L. 



t9 = 1.0, n = 20 



L 


Ref. Price 


RBB price 


RBB var. 


QEP price 


105 


0.029 


0.029 


0.067 


0.029 


106 


00.06 


00.06 


0.165 


00.06 


107 


00.11 


00.11 


0.165 


00.11 


110 


00.43 


00.43 


1.966 


00.43 


111 


00.61 


00.61 


3.022 


00.61 


112 


00.83 


00.84 


4.423 


00.84 


115 


01.77 


01.77 


10.89 


01.78 


120 


04.16 


04.20 


27.81 


04.20 


125 


07.11 


07.14 


45.41 


07.17 


130 


09.87 


09.90 


58.17 


09.92 



Table 5: Up-and-out Call prices from RBB and QEP methods in the local volatility model for r ~ 0.15, 
S = 0.5, 7? = 0.1, cZat 966, T = 1, K = 100, Xq = 100 and for varying values of the barrier L. 
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ODE (15.81 ). there is no hope to improve of price approximations for quantization method by increasing 
the number 7i of discretization steps, so that the RBB method will become more competitive when 
increasing n. To improve the estimations for quantization method we must increase the size of d^r. 
But, this will increase the computation time and the method will become less faster than the RBB one. 
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